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We examine the effect of the trapping potential on the critical temperature, Tc, for the BCS 
transition to a superfluid state in trapped atomic gases of fermions. Tc for an arbitrary power law 
trap is calculated in the Thomas-Fermi approximation. For anharmonic traps, Tc can be increased 
by several orders of magnitude in comparison to a harmonic trap. Our theoretical results indicate 
that, in practice, one could manipulate the critical temperature for the BCS phase transition by 
shaping the traps confining the atomic Fermi gases. 
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Since the observation of Bose-Einstein condensation in 
trapped atomic gases [Q , there has been increasing inter- 
est in the possibility of observing the Bardeen-Cooper- 
Schricffer (BCS) phase transition to the superfluid state 
in dilute fermionic akali gases ||, ^, |5[ ^, 0, ||, ^ |l^, 
|f l| . Currently, experimental efforts in cooling fermionic 
atoms of ^Li [|f| [^] and ^^K @, |l| to the quantum de- 
generate regime have made significant progress, reaching 
temperatures as low as Q.2Tp where Tp is the Fermi tem- 
perature. However, these temperatures are still far be- 
yond the critical temperature required for the BCS phase 
transition, which is at least an order of magnitude lower 

SiilJ. 

On the other hand, the current techniques used to cool 
fermionic atoms lose their utility when the temperature 
becomes much less than Tp. For example, evaporative 
cooling, which has been used for ^"K, becomes ineffec- 
tive for temperatures much less than Tp due to Pauli 
blocking ||l6|] while sympathetic cooling of the fermions 
using bosons, as done with ^Li loses its effectiveness 
when the heat capacity of the bosons falls below that of 
the fermions. Due to these technical obstacles in cool- 
ing fermionic atoms, it seems that one must seek other 
avenues to reach the BCS superfluid state. 

Recently there have been a number of proposals to 
study the possibility of achieving a higher Tc by increas- 
ing the strength of the attractive interactions, i.e., the 
scattering length between the fermions using either a Fes- 
hbach resonance |l^ or photoassociation ||ri|. Since 
Tc is a function of the attractive interaction between 
fermions and the density of states at the Fermi energy, 
Pp, higher transition temperatures can also be achieved 
by increasing the density of states for the fermions at the 
Fermi surface, fn bulk superconductors, the only way to 
modify the density of states is by changing the spatial di- 
mensionality of the system. However, in atomic systems 
the density of states can be controlled through the shape 
and strength of the external trapping potential. 

The ability of an external trapping potential to change 
the critical temperature for a phase transition by mod- 
ifying the density of states, p(e) (see Eq. (^ below), 
was first considered by Bagnato et al. for Bose-Einstein 
condensation . They showed that the critical temper- 



ature for Bose-Einstein condensation, TbeCj ^ ^ func- 
tion of the number of atoms confined in a power law 
trap, Nb, had the power law form Tbec ^ where 
2/9 < rj < 2/3. fn a similar manner, we will show that 
for a trapped gas of Np fermions the density of states at 
the Fermi surface has the form pp ^ Np where [3 varies 
from 1/9 to f/3 depending on the type of trap. Since 
Tc depends exponentially on pp, small changes in (3 can 
have a dramatic effect on Tc. fn all the experiments with 
trapped fermions, the trapping potential is well approxi- 
mated as being harmonic. Our results indicate, however, 
that higher values of Tc/Tp can be obtained using an- 
harmonic power law traps. This indicates that the use 
of anharmonic traps could place Tc within the range of 
current experiments. 

At the ultracold temperatures achieved in fermion ex- 
periments, p-wave collisions between atoms are highly 
suppressed and s-wave collisions between atoms in the 
same internal state are forbidden by the Pauli exclusion 
principle. Since the BCS transition requires an attractive 
interaction between fermions in order to form Cooper 
pairs, the most likely candidate for Cooper pairing in 
alkali gases is an attractive s-wave interaction between 
atoms in different hyperfine states. Fortunately, both 
^Li and appear to be very promising candidates. ^Li 
possesses an anomalously large and negative s-wave scat- 
tering length equal to — 2f60ao |^, where qq is the Bohr 
radius, while for '^'^K, a Feshbach resonance exists for two 
of the hyperfine states which can be used to create the 
required large attractive interaction 

Therefore, we discuss the case of s-wave pairing be- 
tween atoms of mass m in two hyperfine states with 
an s-wave scattering length a < for collisions between 
atoms in different states. Furthermore, we assume that 
the number of atoms in each of the hyperfine states, Np, 
is the same since this is the optimal condition for Cooper 
pairing. In this case, the transition temperature for a 
dilute homogenous gas is given by 

p^ = aexp( ^-V (I) 

Tp \ gPFj 

where g = A'KlT?\a\/m is the coupling constant for the 
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two-body interactions and a = ^/^/tt Ri 0.28 with 

C being the Euler's constant ||, The density of 

states at the Fermi surface, in terms of the Fermi mo- 
mentum hkp, is 

Pp — rnkp /^Ti'^h^ . 

From Eq. we can see that by raising the density of 
states at the Fermi surface (or equivalently, by raising 
kp), there wiU be an exponential increase in the tran- 
sition temperature. Note that even though we restrict 
ourselves to the case of s-wave pairing, our results can 
easily be extended to pairing via higher partial waves 
since in those cases the critical temperature has a similar 



exponential dependence on pp |19 



The goal of this paper is to determine the critical 
temperature for the inhomogeneous Fermi gas where the 
atoms in both hyperfine states are subject to the trapping 
potential. 
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and p, I, and q are positive integers. We refer to the mag- 
nitude of the exponents p, /, and q as the confining power 
of the trap. For simplicity, we assume that the trapping 
potential is independent of the hyperfine state, hence the 
density of fermions is the same for both states, n(r). If 
n(r) varies sufficiently slowly, then we can make a local 
density (or Thomas-Fermi) approximation by assuming 
that at each point in space the Fermi gas can be treated 
as being homogenous ||2^. In that case, a local value 
of the critical temperature, Tb(r), can be calculated us- 
ing Eq.(|]). For temperatures T <^Tp, the local density 
approximation is valid provided the average distance be- 
tween atoms, ~ kp^, is much less than the distances over 
which n(r) changes significantly. Since the changes in 
n(r) are determined by U{r) for which the characteris- 
tic length scales are a, b, and c, we have the condition 
kp^ <C {abcy/^. In addition, for T ^ Tc the coherence 
length for the Cooper pairs, = h^kp /'Kra/S.Q where 
Ao = {■Ke~^)kBTc is the BCS gap at T = 0, should also 
be much less than (abc)^/'^ Q|. 

In the Thomas-Fermi approximation, the chemical po- 
tential, p, is replaced with a local value, /x(r) = /i— J7(r) + 
.gn(r) so that the fermions can be treated at each r as an 
ideal homogenous gas with chemical potential /i(r). The 
term gnir) represents the Hartree-Fock potential expe- 
rienced by each atom due to the atoms in the opposite 
hyperfine state. Since we are interested in temperatures 
T <C Tf, the chemical potential may be approximated 
by the Fermi energy, p ~ Ep = ksTp = h^k %/2 ra, since 
for low temperatures p = Ep + 0{T/Tp)'^ |2l]. Corre- 
spondingly, we can express p{r) in terms of a local Fermi 
wave number, kp{r), defined as 



h^klir) 
2m 



= Ep-U{Y)+gn{v). 
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Using the result n(r) = fc|,(r)/(67r^) for an ideal Fermi 
gas, then gives an expression for the density 
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where A = Gn^ (h^ /2mY' ^ ■ For a dilute gas where 
|ajn(0)^/'^ ^ 1, the Hartree-Fock term in Eqs. (^) and 
(^) can be neglected in a first approximation. The Fermi 
energy is then determined by the requirement that the 
total number of atoms in each spin component be con- 
served. 
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where the integration volume V(Ep) is the volume avail- 
able to a classical particle with total energy Ep, i.e. 
Ep > U{r). Equation (||) can easily be solved for Ep 
for the case of the power law trap (|2|) and gives, 



Ep 



r(i + i/p)r(i + i/or(i + i/g)r(5/2) \8abc 
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where T{x) is the Gamma function and 

S = 3/2+l/p+l/l + l/q. 

Note that the limit p,l,q oo corresponds to a box with 
volume 8abc fl^ . For a harmonic oscillator, p — I — q — 

2, Sabc is equal to 2^^'^lxlytz where li = yjh/mui is the 
harmonic oscillator length along the axis of the trap with 
frequency LOi. Therefore, n = Np/(8abc) can be used to 
define a characteristic atomic density. Note that, with 
the exception of the rigid box, n does not correspond to 
the average density of atoms in the trap. In what follows, 
we will assume that n is fixed and independent of S. 

Equation (^ along with Eq. (S) can be used to calcu- 
late Tc{r). It is clear from Eq. P) that kp{r) is a max- 
imum at the center of the trap and hence, the density of 
states at the local Fermi surface, pp{r) = mkp{Y) /2ir^'h? , 
is a maximum there. Consequently, Tc{y) is largest at 
r = 0. Physically, this means that as the temperature is 
lowered Cooper pairing first occurs at the center of the 
trap and then spreads to the edges of the trap as the 
temperature is lowered still further. Therefore at r = 
and neglecting the Hartree-Fock term in Eq. (||), the 
transition temperature is simply. 
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Evaluating the critical temperature at the center of the 
trap has the added benefit that the Thomas-Fermi ap- 
proximation is expected to be most accurate here. 

From Eqs. (|) and (0) we see that for a given value 
of n, Tc{0)/Tp is an increasing function of 1/5. Fur- 
thermore, Tp = Ep/ks is also an increasing function 
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of 1/5. Therefore, increasing the confining power of the 
trap increases not only the ratio of the critical tempera- 
ture to the Fermi temperature but also the absolute value 
of the critical temperature. More importantly, from the 
perspective of current experiments, anharmonic traps for 
which p,l,q > 5 will result in higher values of Tc (assum- 
ing that all the other terms in Eq. (|^) are the same) . This 
is the central result of this paper. 

To illustrate the effect of varying the confining power of 
the trap, we consider the case oi p — I — q for the values 
of 1 through 5 and oo. We choose values of £j and a, b, 
and c that correspond to parameters used in current ex- 
periments with harmonic traps. In Tables | and || we cal- 
culate the Fermi temperature, Tp, as well as Tc(0) for ^Li 
and *°K, respectively. For ^Li we consider an isotropic 
trap with Np = 10^ and Si ^ hw = h{2TT x lOOs^^), 
which gives a = b — c — ^J^Kjmu} — 5.8/im and 
fi = 6.4 X VS^^cmT^. For ^°K, we use values from 
the experiment by DeMarco and Jin [^. This gives 
Np = 3.5 X 10^ and £i = £2 = = S(27r x 1275"^) and 
£3 = fiiOz — h{2TT X 19. 5s^^). This gives values for the 
characteristic lengths of a = 6 = 2^m and c — 5.09/im 
and a characteristic density of fi — 2.15 x lO^^cm^'^. 
Bohn has recently showed that there exists an exper- 
imentally accessible Feshbach resonance for two of the 
hyperfine states of '^''K that could be accessed to create 
a scattering length of a = — lOOOao 1^. We therefore 
adopt this value for a since in the absence of a Feshbach 
resonance, the scattering lengths for would result in 
unreasonably small values of Tq ■ 

Tables | and Q illustrate the dramatic effect that p has 
on Tp and TcCO). For both ^Li and '"^K, there is a two 
order of magnitude increase in the Fermi temperature as 
p is increased from 1 to cxd. Furthermore, in going from a 
harmonic trap (j) = 2) to a rigid box (p = 00) the Fermi 
temperatures increase by factors of 24 and 36 for ^Li and 
'^''K, respectively. Even more striking is the change in 
Tc{Q) /Tp, which increases by three orders of magnitude 
over the full range of p. Altogether, this implies that by 
increasing the confining power of the trap (p) one could, 
in principle, increase Tc(0) by as much as three orders of 
magnitude in comparison to a harmonic potential. Note 
that in our calculation, we take Ep as the Fermi energy 
of an ideal Fermi gas. Including the attractive atom- 
atom interaction will decrease Ep. However, one must 
now include the Hartree-Fock term, gn{{)), in Eq. (|^), 
which raises the critical temperature. We checked this 
effect numerically and find that for the parameters used 
in the paper, including the atom-atom interactions in- 
creases T'c(O) by a factor of 2 to 3. 

Physically, the increase in Tp and ^^(O) with increas- 
ing confining power is a result of the trap being able 
to confine the atoms to a smaller total volume, which 
thereby increases the local density. For a rigid box, the 
atoms are confined to a volume of 8a6c regardless of 
the value of Ep. For a harmonic trap, the total vol- 
ume occupied by the gas corresponds to the extent of 
the wave function for an atom in the highest occupied 



state with energy Ep ^ f\w, which has a volume much 
larger than the volume of the ground state wave function 
given by ^ abc. In general, the total volume is given by 
V{Ep). For the isotropic form of the power law potential, 
U{v) = £(r/a)P, the volume scales as V{Ep) ~ 7V^/'p+^\ 
This explains the increase in Tp and Tq since, in the limit 
of a homogenous gas, they depend only on the density of 
the fermions. 

Alternatively, the increase in Tc can also be explained 
by examining the density of states for an atom with total 
energy between e and e + de [0 , 




Increasing the confining power of the trap reduces the 
volume of phase space available to an atom with en- 
ergy e, V{t), and as a result, the number of states is 
reduced. In fact, it is easy to show that p(e) '-^ e*~^. 
The reduction in the density of states causes Ep to be 
increased since the fermions are forced to occupy higher 
energy states in order to accommodate all Np atoms. In 
the Thomas-Fermi approximation, the density of states 
on the local Fermi surface at the center of the trap is 
pp{Q) = mkp{0)/2Tr^h'^ ^ \/Ep. Consequently, an in- 
crease in the Fermi energy, Ep, increases the local den- 
sity of states and therefore T'c(O), as can be seen from 
Eq. (0). 

Both Eq. (|l|) and Eq. (g) are valid in the limit of 
a dilute gas, |a|n(0) <C 1. This approximation begins 
to break down for p > h when |a|n(0) > 0.2 for both 
cases considered here. There are, however, two reasons 
why this need not be of any great concern. First, the 
neglect of the Hartree-Fock mean-field in Eq. (||) under- 
estimates the actual density of the gas at the center of 
the trap because the interactions are attractive. The fact 
that attractive interactions increase the local density for 
trapped gases is well known for Bose-Einstein conden- 
sates 1^ and has been previously noted for fermions . 
Secondly, Heiselberg has shown that in the regime of in- 
termediate densities, \a\n > 1, Tc becomes a finite frac- 
tion of Tp with values of Tc/Tp > 0.1 ||. Therefore, 
for large confining powers, our results underestimate the 
actual value of Tc(0). 

Finally we remark that from an experimental point of 
view, changing the confining power of the trapping po- 
tential appears to be realistic. The generation of power 
law traps would most easily be accomplished using far- 
detuned optical dipole traps, which have the added ben- 
efit of producing a confining potential that is indepen- 
dent of the hyperfine state of the atoms. In particular, 
the generation of higher-order Bessel beams with radial 
intensity profiles proportional to Jf{krr), where J/ is a 
Bessel function of integer order and kr is the radial com- 
ponent of the wavevector, have recently been produced 
using an axicon |2^. Bessel beams with / = 1 to 4 and 
radii for the hollow core of the beam on the order of tens 
of micrometers were created. For blue-detuned beams 
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TABLE L Tf and 70(0) as functions of the confining power of tlie isotropic power-law trap, for ''Li. 



p^l = q Ef{1) Tf{K) Tc{0)/Tf rc(0)(K) 

1 1.90 X 10"=*" 1.38 X 10"^ 1.62 x lO"* 2.24 x 10"^ 

2 5.56 X 10-^° 4.03 x 10"^ 3.58 x 10"^ 1.44 x 10" 

3 1.06 X 10"^'^ 7.69 x 10"'' 0.0115 8.84 x 10" 

4 1.63 X 10"^'' 1.18 X 10"'^ 0.022 2.60 x 10" 

5 2.20 X 10"^" 1.59 X 10"® 0.031 4.93 x 10" 
00 1.36 X 10"^* 9.82 x 10"® 0.116 1.14 x 10"' 



TABLE II: Same as Table § for '^"K. 

p^l = q EFjJ) Tf{K) Tc{0)/Tf rc(0)(K) 

1 1.71 X 10"^" 1.24 X 10"'^ 3.93 x 10"* 4.87 x 10"" 

2 5.76 X 10"^° 4.17 x 10"'^ 7.8 x 10"^ 3.26 x 10"^ 

3 1.19 X 10"^^ 8.62 X lO"'^ 0.023 1.98 x 10"* 

4 1.93 X 10"^'' 1.40 X 10"® 0.04 5.54 x 10"* 

5 2.72 X 10"^^ 1.97 x 10"® 0.054 1.06 x 10"^ 
00 2.12 X 10"^* 1.54 X lO"-'^ 0.155 2.37 x 10"® 



and small fc^r, the radial optical dipole potential expe- 
rienced by the atoms is proportional to (fc^r)^'. Two 
perpendicular intersecting Bessel beams with I > 1 could 
be used to create an anharmonic potential. It is worth 
noting that evaporative cooling of a two-component gas 
of ^Li to temperatures below Tp in an optical trap has 
recently been demonstrated experimentally . 

In conclusion we have examined the effect that a 
power-law trapping potential has on the BCS transition 
temperature. We have shown that by increasing the con- 



fining power of the trap, one can obtain values of Tc 
that are several orders of magnitude larger than the cor- 
responding harmonic trap. The origin of the increase 
is the ability of tighter traps to confine the atoms to a 
smaller total volume. 
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